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We study diffuse gamma-ray emission at intermediate Galactic latitudes measured by the Fermi Large Area 
Telescope with the aim of searching for a signal from dark matter annihilation or decay. In the absence of a 
robust dark matter signal, constraints are presented. We set both, conservative dark matter limits requiring that 
the dark matter signal does not exceed the observed diffuse gamma-ray emission and limits derived based on 
modeling the foreground astrophysical diffuse emission. Uncertainties in several parameters which characterize 
conventional astrophysical emission are taken into account using a profile likelihood formalism. The resulting 
limits impact the range of particle masses over which dark matter thermal production in the early Universe is 
possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of 
dark matter. 



1. Introduction 



Most of the mass in our Universe is in the form of 
yet un-identified particles (i.e. Dark Matter (DM)) 
which have been detected only through their gravita- 
tional interactions thus far. In one of the most at- 
tractive frameworks to explain the DM problem (the 
WIMP paradigm) those particles are expected to self 
annihilate to stable standard model particles, produc- 
ing gamma rays, electrons and protons. Due to our 
proximity to the center of the Milky Way DM halo, 
such gamma ray emission originating in our Galaxy 
would appear as a diffuse signal. 

At the same time, the majority of the Galactic dif- 
fuse emission is produced through radiative losses of 
cosmic-ray (CR) electrons and nucleons in the inter- 
stellar medium. Modeling of this emission presents 
one of the major challenges when looking for subdom- 
inant signals from dark matter. 

In this analysis we test the diffuse LAT data for a 
contribution from the DM signal by performing a fit of 
the spectral and spatial distributions of the expected 
photons at intermediate Galactic latitudes. In doing 
so, we take into account the most up-to-date modeling 
of the established astrophysical signal, [Ackcrma nn ei\ 



at] [2012|b 



Our aim is to constrain the DM prop- 
erties and treat the parameters of the astrophysical 
diffuse gamma-ray background as nuisance parame- 
ters. Those parameters are typically correlated with 
the assumed DM content and it is thus important 
to scan over them together with the DM parameter 
space, since they affect directly the DM fit. 

Besides this approach, we will also quote conserva- 
tive upper limits using the data only (i.e. without 
performing any modeling of the astrophysical back- 
ground) . 



2. Modeling of the high-energy Galactic 
diffuse emission 



Wc follow Ackcrmann et al. 



GALPROP code v54 Strong et al 



2012b 



2000 



in using the 
to calculate 

the propagation and distribution of CRs in the Galaxy 
and the whole sky diffuse emission, as well as the sig- 
nal from DM. Several parameters enter the CR prop- 
agation modeling, see Ackermann et al. [2012b| for 
more detail: 



the 



half- 



the distribution of CR sources, 
height of the diffusive halo Zh, the radial extent of 
the halo Rh, the nucleon and electron injection spec- 
trum, the normalization of the diffusion coefficient 
D , the rigidity dependence of the diffusion coeffi- 
cient 5, (D(p) = D (p/p )~ 5 with p being the refer- 
ence rigidity) and the Alfven speed va, (parametrizing 
the strength of re-acceleration of CRs in the ISM via 
Alfven waves) and the velocity of the Galactic winds 
perpendicular to the Galactic Plane V c . Interactions 
of the CRs with the interstellar medium (ISM) and 
interstellar radiation field (ISRF) produce three dis- 
tinct components of the gamma-ray emission: photons 
from the decay of neutral pions produced in the inter- 
action of the CR nucleons with the interstellar gas, 
bremsstrahlung of the CR electron population on the 
interstellar gas and their inverse Compton scattering 
off the interstellar radiation field. 

various standard pa- 



in Ackermann et al. 2012b 



rameters of the CR propagation were studied in a fit 
to CR data and it was shown that they represent well 
the gamma-ray sky, although various residuals (at a 



30% level Ackermann et al. 2012b ), both at small 



and large scales, remain. These residuals can be as- 
cribed to various limitations of the models: imperfec- 
tions in the modeling of gas and ISRF components, 
simplified assumptions in the propagation set-up, un- 
resolved point sources, and large scale structures like 
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Loop I Casandjian et al. 2009 or the Galactic Bub- 
bles Su et al. |2010| . Since residuals do not seem ob- 
viously related to DM, we focus in the following on 
setting limits on the possible DM signal, rather than 
searching for a DM signal. 

In our work, we use the results of the fits to the CR 
data from Ackermann et al. 2012b but we allow for 



more freedom in certain parameters governing the CR 
distribution and known astrophysical diffuse emission 
and constrain these parameters by fitting the models 
to the LAT gamma-ray data. 



3. DM maps 

Numerical simulations of Milky Way size halos re- 
veal a smooth halo which contains large number of 



subhalos Diemand et al. 2007 , |Springel et al. 2008 



The properties of the smooth halo seem to be well 
understood, at least on the scales resolved by simu- 
lations, while the properties of the subhalo popula- 
tion are more model dependent. In the inner < 20° 
region of the Galaxy, the smooth component is ex- 



pected to dominate, Diemand et al. 2007b , Springel 



et al. 2008b , Pieri et al. 2011 and we conservatively 



consider only the smooth component in this work. 
We parametrize the smooth DM density p with a 



NFW spatial profile Navarro et al. 1996] 



Pi> 



PqR s 

r(l + r/R s f 



(1) 



and a cored (isothermal-sphere) profile Bcgcman et al. 

ImT 



p(r) 



Po (Rq + R 2 C 

(r 2 + Rl) 



For the local densi ty of DM we take the value of 
p = 0.43 GeV cm" ~ 



Salucci et al. 2010 , and the scale 
radius of R s = 20 kpc (for NFW) and R c = 2.8 kpc 
(isothermal profile). We also set the distance of the 
solar system from the center of the Galaxy to the value 
Rq = 8.5 kpc. For the annihilation/decay spectra we 
consider three channels with distinctly different sig- 
natures: annihilation/decay into the bb channel, into 
p + \x~ , and into t + t~. In the first case gamma rays 
are produced through hadronization of annihilation 
products and subsequent pion decay. The resulting 
spectra are similar for all channels in which DM pro- 
duces heavy quarks and gauge bosons and this chan- 
nel is therefore representative for a large set of particle 
physics models. The choice of leptonic channels pro- 
vided by the second and third scenarios, is motivated 



by the dark matter interpretation Grasso et al. 



2008 



of the PAMELA positron fraction Adriani et al 
and the Fermi LAT electrons plus positrons | A 
al. 2009 measurements. In this case, gamma rays are 



2009 



ido et 



dominantly produced through radiative processes of 
electrons, as well as through the Final State Radiation 
(FSR). We produce the DM maps with a version of 
GALPR0P slightly modified to implement custom DM 
profiles and injection spectra (which are calculated by 
using the PPPC4DMID tool described in |Cirelli et ~aT\ 
2010 and include a contribution from electro-weak 



bremsstrahhmg) . 



4. Approach to set DM limits 

We use 24 months of LAT data in the energy range 
between 1 and 100 GeV (but, we use energies up to 
400 GeV when deriving DM limits with no assumption 
on the astrophysical background) . We use only events 
classified as gamma rays in the P7CLEAN event se- 
lection and the corresponding P7CLEAN_V6 instrument 
response functions (IRFs) 1 . Structures like Loop I and 
the Galactic Bubbles appear mainly at high Galac- 
tic latitudes and to limit their effects on the fitting 
we will consider a ROI in Galactic latitude, 6, of 
5° < \b\ < 15°, and Galactic longitude, I, \l\ < 80°. 
We mask the region \b\ < 5° along the Galactic Plane, 
in order to reduce the uncertainty due to the modeling 
of the astrophysical and DM emission profiles. 

4.1. DM limits with no assumption on the 
astrophysical background 

To set these type of limits we first convolve a given 
DM model with the Fermi LAT instrument response 
functions (IRFs) to obtain the counts expected from 
DM annihilation. The expected counts are then com- 
pared with the observed counts in our ROI and the 
upper limit is set to the minimum DM normalization 
which gives counts in excess of the observed ones in 
at least one bin, i.e. we set 3cr upper limits given by 
the requirement uidm — 3\/niDM > n ij where uhjm 
is the expected number of counts from DM in the bin 
i and n, the actual observed number of counts. 



4.2. DM limits with modeling of 
astrophysical background 

In this analysis we model the diffuse emission as 
a combination of a dark matter and a parameterized 
conventional astrophysical signal and we derive the 
limits on the DM contribution using the profile likeli- 
hood method. 

The parameters we use to describe conventional as- 
trophysical emission and their ranges are summarized 
in Table HI 



1 http://fermi.gsfc. nasa.gov/ssc/ 
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For each DM channel and mass, the model which 
describes the LAT data best is determined by maxi- 
mizing the likelihood function defined as 



Ld 



'DM) 



Jk{VDM,Ct) 



max s Y\_Pik{ni;d, 9 



DM) 



(3) 

where i are spatial and spectral bins, P^ is the Pois- 
son distribution for observing n, events in bin i given 
an expectation value that depends on the parameter 
set (6dm, cT). 9dm is the intensity of the DM com- 
ponent, a represents the set of parameters which en- 
ter the astrophysical diffuse emission model as linear 
pre-factors to the individual model components (cf. 
equation [4] below), while k denotes the set of parame- 
ters which enter in a non-linear way. We sample non- 
linear parameters on a grid in the k parameter space 
(see Table III]) . On each point on a grid of non-linear 
parameters and for each fixed value of DM normal- 
ization 6dm we hud the values of linear parameters 
(a) which maximize the likelihood and the value of 
the likelihood itself at the maximum, thus determin- 
ing the curve L^{0dm)- The outlined procedure is 
then repeated for each set of values of the non-linear 
propagation and injection parameters to obtain the 
full set of profile likelihood curves. We scan over the 
three parameters: electron injection index, the height 
of the diffusive halo and the gas to dust ratio which 
parametrizes different gas column densities (see Table 
|llj ). In this way we end up with a set of k profiles of 
likelihood £fc(#DJVf)j one f° r each combination of the 
non-linear parameters. The envelope of these curves 
then approximates the final profile likelihood curve, 
L{0dm)i where all the parameters, linear and non- 
linear have been included in the profile. Limits are cal- 
culated from the profile likelihood function by finding 
the DM ,lim values for which L(9 DM ,Um)/L(9DM,max) 
is exp(— 9/2) and exp(— 25/2), for 3 and 5 a C.L. lim- 
its, respectively. 

In each step, the maps (produced by the GALPROP 
code) which are used for fits, and their normalization 
parameters are combined as: 

i *■ j 

j 
a x (x 7 + X* c ) + J] a iGB, m IGB m . (4) 



here, the sum over i is the sum over all step-like 
distributions of cosmic ray source (CRSD) functions 2 , 



W V brcmss 



the sum over j corresponds to the sum over all Galac- 
tocentric annuli (details of the procedure of a place- 
ment of the gas in Galactocentric annuli and their 
boundaries are given in jAckermann et~aT. 2012b ). H 
denotes the gamma-ray emission from atomic and ion- 
ized interstellar gas while H2 the one from molecular 
hydrogen 3 and IC the Inverse Compton emission. \ 1 
and Xic are the prompt and Inverse Compton (when 
present) DM contribution and a x the overall DM nor- 
malization. Additionally an isotropic component aris- 
ing from the extragalactic gamma-ray background and 
misclassified charged particles needs to be included to 
fit the Fermi LAT data. aiGB,m and IGB m denote 
the Isotropic Gamma-ray Background (IGB) intensity 
for each of the five energy bins over which the index 
?n runs. In all the rest of the expression the energy 
index m is implicit since we don't allow for the free- 
dom of varying the GALPROP output from energy bin 
to energy bin. We do not include sources in the fit as 
we use a mask to filter point sources from the 1FGL 



catalog Abdo et al. 2010 



5. Results 

An important point to note about our fitting proce- 
dure is that, for each DM model, the global minimum 
we found lies within the 3(5) a regions of many dif- 
ferent models, providing a check against a bias in our 
procedure. This point is illustrated in Figurefll where 
the profile likelihoods for the three nonlinear param- 
eters, Zh, 7e,2 and d2HI, are shown. To ease read- 
ing of the figure the profiling is actually performed 
with further grouping DM models with different DM 
masses, but keeping the different DM channels, DM 
profiles and the annihilation/decay cases separately. 
The curve for the fit without DM is also shown for 
comparison. Each resulting curve has been further 
rescaled to a common minimum, since we are in- 
terested in showing that several models are within 
— 2AlogZ/ < 25 around the minimum for each DM fit. 



2 CRSD over the Galactic plane is a critical parameter for 



DM searches. CRSDs are traditionally modeled from the direct 
observation of tracers of SNR and can be observationally biased. 
In order to circumvent this problem, we define a parametric 
CRSD as sum of step functions in Galactocentric radius R, and 
treat the normalization of each step as a free parameter in the 
fit to gamma rays. In addition, in order to have conservative 
and robust limits we set to zero the e,p CRSDs in the inner 
Galaxy region, within 3 kpc of the Galactic Center. In this 
way, potential e and p CR sources which would be required 
in the inner Galaxy will be potentially compensated by DM, 
producing conservative constraints. 

3 It should be noted that in our case, where we mask ±5° 
along the plane, the expression actually simplifies considerably 
since only the local ring Xqq factor enters the sum, since all 
the other H2 rings do not extend further than 5 degrees from 
the plane. 
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The 7e,2 profile, for example, indicates that all mod- 
els with 7 e: 2 from 1.9 to 2.4 are within — 2AlogL < 25 
around the minimum illustrating that the sampling 
around each of the minima for the six DM models is 
dense. The Zh profile indicates that basically all the 
considered values of Zh are close to the absolute min- 
ima. This last result is not surprising since, within our 
low-latitude ROI, we have little sensitivity to different 
Zh and basically all of them fit equally well. There is 
some tendency to favor higher values of Zh when DM 
is not included in the fit, while with DM the trend is 
inverted. Although the feature is not extremely sig- 
nificant it is potentially very interesting. 

Upper limits on the velocity averaged annihilation 
cross section into various channels are shown in Fig.[2j 
for isothermal profile of the DM halo 4 , together with 
regions of parameter space which provide a good fit to 
PAMELA (purple) and Fermi LAT (blue) CR electron 



and positron data Circlli et al. 2009 



The resulting DM limits are comparable with the 
limits from LAT searches for a si gnal from DM an- 
nihilation/decay in dwarf galaxies Ackermann et al. 
2011 . In particular, as shown in Figure [2J for masses 



around 20 GeV the thermal relic value of the anni- 
hilation cross section is reached, both for the bb and 
t + t~ channels. 

In addition to the parameters listed in Table |LT| we 
check the importance of uncertainties in additional as- 
trophysical parameters, but in a more simplified setup: 
we set a particular model as reference and then we 
vary each parameter one at a time, keeping the others 
fixed, and for each case we calculate the percentage 
variation in DM limits for selected DM models. We 
find that the produces chan g es in the limits of less 

Overall, rather 



than 10% Ackermann et al. 2012 



than being due to residual astrophysical model uncer- 
tainties, the remaining major uncertainties in the DM 
constraints from the Halo region come from the mod- 
eling of the DM signal itself. The main uncertainty is 
in the normalization of the DM profile, which is fixed 
through the local value of the DM density. We use 
the recent determination p = 0.43 GeV cm~ 3 from 
ilucci et al. 2010 , which has, however, a large uncer- 
tainty, with values in the range 0.2-0.7 GeV cm -3 still 
viable. A large uncertainty in p$ is particular impor- 
tant for annihilation constraints since they scale like 
p\, while for constraints on decaying DM the scaling 
is only linear. A less important role is played by the 
uncertainties in the DM profile, since in our region of 
interest different profiles predict similar DM densities. 



References 

M. Ackermann et al. [The Fermi-LAT collaboration], 
Astrophys. J. 761, 91, 2012. 

M. Ackermann et al. [The Fermi-LAT collaboration], 
Astrophys. J., 750, 3 (2012), 

A. W. Strong, I. V. Moskalenko and O. Reimer, 
Astrophys. J. 537, 763 (2000), url: 
http://galprop. Stanford, edu/webrun.php. 

J. -M. Casandjian, I. Grenier [The Fermi-LAT collab- 
oration], [arXiv:0912.3478| [astro-ph.HE]]. (2009) 

M. Su, T. R. Slatyer, D. P. Finkbeiner, Astrophys. J. 
724, 1044-1082 (2010). 

J. Diemand, M. Kuhlen and P. Madau, Astrophys. J. 
667 (2007) 859 

V. Springel et al., Mon. Not. Roy. Astron. Soc. 391, 
1685 (2008) 

J. Diemand, M. Kuhlen and P. Madau, Astrophys. J. 
657 (2007) 262 

V. Springel, S. D. M. White, C. S. Frenk, 
J. F. Navarro, A. Jenkins, M. Vogelsberger, J. Wang 
and A. Ludlow et al, |arXiv:08 09.0894 [astro-ph] 
(2008). 

L. Pieri, J. Lavalle, G. Bertone and E. Branchini, 
Phys. Rev. D 83, 023518 (2011) 

J. F. Navarro, C. S. Frenk and S. D. M. White, As- 
trophys. J. 462, 563 (1996). 

K. G. Begeman, A. H. Broeils and R. H. Sanders, 
Mon. Not. Roy. Astron. Soc. 249, 523 (1991). 

P. Salucci, F. Nesti, G. Gentile, C. F. Martins, Astron. 
Astrophys. 523, A83 (2010). 

D. Grasso et al. [FERMI-LAT Collaboration], 
Astropart. Phys. 32, 140 (2009). P. Meade, 
M. Papucci, A. Strumia and T. Volansky, Nucl. 
Phys. B 831, 178 (2010). L. Bergstrom, J. Edsjo 
and G. Zaharijas, Phys. Rev. Lett. 103, 031103 
(2009). 

O. Adriani et al. [PAMELA Collaboration], Nature 
458, 607 (2009). 

A. A. Abdo et al. [The Fermi LAT Collaboration], 
Phys. Rev. Lett. 102, 181101 (2009). 

M. Cirelli et al., JCAP 1103 (2011) 051. 

A. A. Abdo et al. [The Fermi-LAT collaboration], As- 
trophys. J. Suppl. 188 (2010) 405-436. 

M. Circlli, P. Panci, P. D. Serpico, Nucl. Phys. B840 
(2010) 284-303. 

M. Ackermann et al. [Fermi-LAT Collaboration], 
Phys. Rev. Lett. 107, 241302 (2011). 



4 Limits obtained using the NFW profile are only slightly 
better. 
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Table I Paramateres 



Non linear Parameters 


Symbol 


Grid values 


index of the injection CRE spectrum 

half height of the diffusive halo 

dust to HI ratio 


7e,2 
Zh 

d2HI 


1.800, 1.925, 2.050, 2.175, 2.300, 2.425, 2.550, 2.675 

2, 4, 6, 8, 10, 15 kpc 

(0.0120, 0.0130, 0.0140, 0.0150, 0.0160, 0.0170) xlO -20 mag cm 2 


Linear Parameters 


Symbol 


Range of variation 


eCRSD and pCRSD coefficients 

local Ebto CO factor 

IGB normalization in various energy bins 

DM normalization 


\rloc 

A co 

CtIGB,m 

a x 


0,+oo 

0-30 xlO 20 cm" 2 (K km s" 1 )" 1 

free 

free 



Table II Summary table of the parameters varied in the fit. The top part of the table show the non linear parameters 
and the grid values at which the likelihood is computed. The bottom part show the linear parameters and the range of 
variation allowed in the fit. The coefficients of the CRSDs are forced to be positive, apart c\' v and c^ p which are set to 
zero. The local Xco ratio is restricted to vary in the range 0-30 xlO 20 cm" 2 (K km s -1 ) -1 , while aiGB,m and a x are 
left free to assume both positive and negative values. See the text for more details. 





0.014 0.015 

cteHltur^magcm 2 ; 



bb. DEC, NFW 
mu, DEC, NFW 
tau, DEC, NFW 




z [kpc] 



Figure 1: Profile likelihood curves for Zh, 7e,2 and d2HI. The various curves refer to the case of no DM or different DM 
models (see the legend in the figure, where we mark a dominant decay (DEC) or annihilation (AN) channel and the 
assumed DM profile). All minima are normalized to the same level. Horizontal dotted lines indicate a difference in 
— 2AlogL from the minimum of 9 (3<r) and 25 (5a). 
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■ FSR, w/o background modeling 

1 IC+FSR, constrained free source fits 



30- 
5a- 



10 



10 2 10 3 10 4 10 10 2 10 3 10 4 

Figure 2: Upper limits on the velocity averaged DM annihilation cross-section (left) and decay time (right) including a 
model of the astrophysical background compared with the limits obtained with no modeling of the background. Limits 
are shown for bb (upper), n + fi~ (middle) and t + t~ (lower panel) channels, for a DM distribution given by the 
isothermal distribution. The horizontal line marks the thermal decoupling cross section expected for a generic WIMP 
candidate. The regions of parameter space which provide a good fit to PAMELA (purple) and Fermi LAT (blue) CR 
electron and positron data are also shown. 
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